Astronomy & Astrophysics manuscript no. 1 1407 


© ESO 2009 


March 2, 2009 





(N 



Received / Accepted 

ABSTRACT 
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$-H ' Context. A short duration burst reminiscent of a soft gamma-ray repeater/anomalous X-ray pulsar behaviour was detected in the 

' direction of LS I -1-61 303 by the Swift satellite. While the association with this well known gamma-ray binary is likely, a different 

I origin cannot be excluded. 

I ' , Aims. We explore the error box of this unexpected flaring event and establish the radio, near-infrared and X-ray sources in our search 

for any peculiar alternative counterpart. 
' Mettiods. We carried out a combined analysis of archive Very Large Array radio data of LS 1 -1-61 303 sensitive to both compact 

^ I and extended emission. We also reanalysed previous near infrared observations with the 3.5 m telescope of the Centro Astronomico 

, Hispano Aleman and X-ray observations with the Chandra satellite. 

Results. Our deep radio maps of the LS 1 +61 303 environment represent a significant advancement on previous work and 16 compact 
radio sources in the LS 1 +61 303 vicinity are detected. For some detections, we also identify near infrared and X-ray counterparts. 
I Extended emission features in the field are also detected and confirmed. The possible connection of some of these sources with the 

, observed flaring event is considered. Based on these data, we are unable to claim a clear association between the Swift-BAT flare 

■ and any of the sources reported here. However, this study represents the most sophisticated attempt to determine possible alternative 

' counterparts other than LS 1 +61 303. 

(<~^ I Conclusions. 

^ I Key words. X-ray: stars - Radio continuum: stars - Infrared: general - X-rays: binaries - Gamma rays: observations 

X: , 

^ 1. Introduction within its 15-150 keV energy range pe Pasquale et al. 2008) 1. 

■ ■ ■ The calculated location of this soft gamma-ray short-duration 

LS I +61 303 (V615 Cas) is a gamma-ray binary origmally burst-event was found to have an eiTor of 2:2 and the po- 

discovered m 1977 in the radio du ring a survey for vari- ^^^^^^ lS I +61 303 was found to be clearly consistent 



able sources in th e Galactic plane ( |Uregory & Taylor 1978[ ^-^^ j^is event (|Barthelmy et al. 2008|. Given this coincidence 



Gregory etal. 1979). The orbital period is about 26.5 d t^is burst was associated with magnetar-like activity linked 



( |Hutchings&Cramptonl981p and it has been detected in data ^ y^^^g highly magnetized pulsar in the binary system 
of frequencies between radio fiaylor & Gregor y 19821 1984) dpubus & Giebels 200 8). Unusual X-ray activity (hard high flux 



and high energy gamma-rays (Albert et al. 2008 ). The physical gpQ^^ j^^^ ^1^^ ^^^^ ^gp^j.^^^ iRay et aL (2008)1 just a 
interpretation of the LS I +61 303 emission across the com- ^^^y.^ ^^^^g ^^^^^ p^A data from the TJXre satellite. If 
plete electromagneti c spectrum stiU remains a matter of debate QpQ^ originate in LS I +6 1 303, an accretion disk would be 
d Romeroetal. 20 03- Two different models have been proposed necessary to explain the nature of this source. However, the PCA 
to explain flie full spectral energy dist ribution: (i) a micro- field of view of ~ 1 degree includes many additional sources that 
quasar X-ray binary (| Bosch-Ramon e IaL2006 ), and (ii) a non- ^^^j^ responsible for the QPO behaviour, 
accreting pulsar interacting with the envelope of the rapidly ro- 
tating Be star dPubus 20061) . Despite these facts, one cannot exclude the possibil- 
On 2008 September 10, the Swift Burst Alert Telescope ity of another unrelated source being responsible for the 
(BAT) detected a burst in the direction of LS I +61 303 observed gamma-ray flare. A population of faint X-ray 

sources in the vicinity of LS I +61 303 were reported by 

Send offprint requests to: A. J. Mufioz-Arjonilla |Rea & Torres (2008)] who suggested that one of these sources 
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Table 1. Radio, near infrared, and X-ray observations of 
LS I +61 303 field used in this paper. 





Instrument 


Date 


Band 


Integ. 
time 






08-Jun-1992 


C (6 cm) 


9.8 h 


Radio 


VLA CnD 


27-Jun-1992 


L (20 cm) 


3.0 h 






09-Sep-1993* 


C (6 cm) 


7.8 h 






13-Sep-1993* 


C (6 cm) 


7.9 h 


Infrared 


CAHA 3.5m 


25-Sep-2007 
25-Sep-2007 


/ (1.25 /^m) 
H (1.65 /im) 


905 s 
905 s 






25-Sep-2007 


K, (2.18 /^m) 


1811 s 


X-rays 


Chandra 


08-Apr-2006 


0.5-7.0 keV 


49.9 ks 



(*) Already used by M98. 



might be the quiescent counterpart of a new transient magne- 
tar. Munoz-Arjonilla et al. (2008)lalso reported additional X-ray 
and radio sources coincident with the 5w//f-BAT error circle and 
some of them with stellar- like counterparts. 

We present extensive radio. X-ray, and near infrared obser- 
vations of the ^wi/f-BAT error circle at the location of this 2008 
September 10 event. Both archival data and observations con- 
ducted by the authors were used in this paper as described in the 
log in Table [1] Different populations of sources were detected 
and their main observational properties are reported. A few in- 
teresting objects in the direction of the magnetar-like flare are 
highlighted and the possibility of them being alternative candi- 
date counterparts is assessed. The census of radio/X-ray sources 
reported here represents the most complete study to date of al- 
ternative counterpart candidates to the 5w/ff-BAT event. 



2. Radio sources within tlie Swift-EK\ error circle 

The magnetar-hke event that renewed interest in studying 
LS I +61 303 is probably related to a compact object formed 
as a result of a supernova event. Therefore, both compact and 
extended radio features could play a role in our understanding of 
this phenomenon. The environment of LS I +61 303 in the radio 
was studied by Ma rtf et al. (1998)| (hereafter M98) at 6 cm wave- 
length using the Very Large Array (VLA) of the National Radio 
Astronomy Observatory (NRAO) with the array in its CnD con- 
figuration providing appropriate sensitivity to both compact and 
extended sources. We developed the M98 approach improving 
their sensitivity by combining additional 6 cm VLA archive data 
acquired by the same compact array configuration (see Table [TJ 
to a total of 25.5 h of on-source time. 

The AIPS package of NRAO was used for the standard in- 
terferometer data processing and self-calibration. We removed 
the variable LS I +61 303 core to avoid artifacts and replaced 
it for cosmetic reasons with a constant point source compo- 
nent with the observed average flux density. We also removed a 
nearby bright radio source, whose presence affected significantly 
the dynamic range of the maps. The final result is presented in 
Fig.[ri where extended emission is enhanced with a slight taper. 
Compact sources are more clearly evident in the non-tapered ra- 
dio map of Fig. |2] with a rms noise of 13 yuJy beam"', signifi- 
cantly lower than that of M98. The observational properties of 
these compact radio sources are listed in Table|2] where J2000.0 
positions are given. 
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Fig. 1. Radio map of the LS I +61 303 environment at the 6 cm 
wavelength obtained after combining all VLA observations dis- 
cussed in the text. The visibility data has been slightly tapered 
with a 15 k/l Gaussian in order to enhance the extended emission 
in the field. The X-ray binary LS I +61 303 and the nearby plane- 
tary nebula G135.6+01.0 (B ond et al. 2003)] are labefled and the 
other relevant sources numbered in right ascension order. Radio 
sources marked with grey numbers are better seen in the non- 
tapered map of Fig. |2] The rms noise is at the level of 14 yuJy 
beam"'. The beam size is shown at the bottom left corner and 
corresponds to 17". x 15". 7, with position angle of -53°. 3. 
Contours shown are -3, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 15, 
20, 25, 50, 100, 200, 500, 1000, 2000 and 3000 times the rms 
noise. 



3. Near infrared observations of tlie ^wZ/i'-BAT error 
circle 

The field of LS I +61 303 was observed at near infrared wave- 
lengths with the 3.5 m telescope at the Centro Astrondmico 
Hispano Aleman (CAHA) in Almeria (Spain), one year be- 
fore the Swift event. The OMEGA2000 camera was used and 
the images were taken through the J, H and Ks filters. CAHA 
observations were processed following the standard proce- 
dures for sky background subtraction, flat-fielding and median- 
combining of individual frames using the IRAFQ software pack- 
age. Astrometry in the final frames was determined by identi- 
fying about twenty stars in the field for which positions were 
retrieved from the 2MASS catalog. The relevant part of the re- 
sulting image is shown in the right panel of Fig.|2] Photometric 
and astrometric data derived from this image are included in 
Table |2] for VLA radio sources with a near infrared counter- 
part. The errors in the photometric measurements are domi- 
nated by the uncertainty in the zero point, which is estimated 
to be about 0. 1 mag in all filters. The existence of a candidate 
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RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Fig. 2. Left. Zoomed contour radio map of the LS I +61 303 environment at 6 cm after combining all the VLA data used in this 
paper Natural weight with no taper has been used. The big circle corresponds to the refined 90% confidence BAT location for the 
magnetar-like flare in the direction of LS I +61 303. The rms noise is at the level of 13 juJy beam The beam size is shown at 
the bottom right corner and corresponds to 12". 7 x 9". 9, with position angle of -58°. 8. Contours shown are -3, 3, 4, 5, 6, 7, 8, 9, 
10, 11, 12, 13, 15, 20, 25, 50, 100, 200, 500, 1000, 2000 and 3000 times the rms noise. Crosses show the X-ray Chandra ACIS-I 
sources with a signal-to-noise ratio above 3. Right. The same field as observed in the near infrared ^T^-band with the CAHA 3.5 m 
telescope and the OMEGA2000 camera. 



counterpart was assessed based on the statistical parameter r 
( |Allington-Smith et al. 1982^ defined to be; 

I (Aa cos 5)2 (Mf 

where Aa and A6 are the differences between the measured near 
infrared and radio positions, and (Ta,iR, cr^rad, cr^.iR and crs,md 
are the infrared and radio uncertainties in right ascension and 
declination. 

The offsets between radio and near infrared positions are also 
listed in Table |2l A value of r < 3 is taken to be indicative of 
astrometric coincidence within the combined radio and near in- 
frared position errors. 

4. X-ray sources within the Swift-BAT error circle 

We also used X-ray data from an observation of the LS I +6 1 303 
field obtained with Chandra about two years before the Swift 
event using the standard ACIS-I setup in VF mode during a total 
of 49.9 ks. The data were reduced using the Chandra Interactive 
Analysis of Observations software package (CIAO v4.0) and 
CALDB V3.4.5. The analysis of LS I +61 3 03 data itself was dis- 
cussed and its results reported in detail by |Paredes et al. (2007)| 
We used the CIAO tool wavdetect to obtain a list of can- 
didate sources with data of a signal-to-noise ratio above 3cr. 
Most of the sources were very faint and the counts insuffi- 
cient for a robust spectral analysis, so we obtained an exposure 



map of the observation and derived the photon fluxes and hard- 
ness ratios of the detected sources. In Table [3] we present the 
fluxes in the 0.5-7.0 keV band and the hardness ratio, defined as 
(H - S)/{H + S), where H is the flux in the 2.0-7.0 keV band 
and S is the flux in the 0.5-2.0 keV band. 

The X-ray observation completely covers the area observed 
in radio with the VLA. In this region, we found a total of 40 
sources, excluding LS I +61 303 itself, 39 of which were previ- 
ously unknown. 



5. Discussion and conclusions 

Our deep radio maps in Figs. [1] and |2] exhibit 16 compact radio 
sources in addition to LS I +61 303 with a peak flux density 
of a factor of four above the rms noise. When we restricted our- 
selves to the Swift-BAJi error circle, three compact radio sources 
(labelled as # 03, 09 and 1 1 in Table |2]l appear to be located 
within or close to the refined BAT location. We measured their 
flux densities at 6 cm at the three epochs quoted in Table[T]sepa- 
rately, and no evidence of variability was found in any source. 
The VLA sources # 03 and 09 are consistent with point-like 
near infrared counterparts (see Table |2]i. Source # 03 coincides 
with a particularly bright object (K,, - 10.0) and its observed 
colours (e. g. J - K, - 1.4) are reminiscent of a giant star with a 
late spectral type, provided that interstellar extinction is not too 
high. Source # 09 appears to be a highly reddened object with 
J - Ks - 2.5 mag and its stellar or extragalactic nature cannot 
be established based on our photometry alone. Neither of them 
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Table 2. Compact radio sources in the field around LS I +61 303 and their candidate near infrared counterparts 



VLA 


Right Ascension 


Declination 


Peak flux dens. 


Total flux dens. 


Aa cos 5 


AO 


r 


J 


H 




Id. # 


(hms) 


(° ' ") 


(mjy beam"') 


(mJy) 


(") 


(") 




(mag) 


(mag) 


(mag) 


QJ*(a) 


02 39 59.44(0.04) 


+61 17 03.3(0.2) 


0.87(0.04) 


1.00(0.07) 


+0.43 


+0.00 


0.71 


20.1 


18.8 


18.1 


02 


02 40 05.63(0.13) 


+61 15 13.5(1.0) 


ATT / A A /I \ 

0.23(0.04) 


A y1 1 /A AA\ 

0.41(0.09) 














03* 


U2 4U lz.2/(U. 11) 


1 1 1/1 oo /I /"A n\ 
+0l 14 22.4(U. / ) 


A OA/A r\A\ 

0.20(0.04) 


A 1 C/A AC\ 

O.Ij(O.Oj) 


1 A 1 

-l.Ul 


+0. /O 


1 IT 
1.1 / 


11/1 
1 1.4 


1 A T 
10.3 


1 A A 
10.0 


04 


0/ 40 12. / /(0.24) 


1 \n OC A/1 A\ 

+61 1 / 25.9(1.4) 


All /A A /I \ 

0.11(0.04) 


A 1 A/ A A^\ 

0.10(0.06) 


-0.94 


-0.10 


0.27 




18.8 


17.7 


n< 
Uj 


AO /1A QO/A /n\ 


+Ol Ij ^i.L\LA) 


A AO/A A/1 \ 

U.US(U.U4) 


A 1 T/A 1 A\ 
U.l /(U.IU) 


1 A ^A 

+U.jU 


A AA 


A Q 
U.JO 


1 A 

lo.U 


1 "7 1 
t /.I 


1 A A 

io.u 


Do 


AO /I A OT Q/1/A 

U2 4U 2/.i4(U.34) 


1 1 1 c 1 o 1 A\ 
+61 15 18.6(1.9) 


A Cil /r\ A/1 \ 

U.U/(U. 04) 


A AiC/A A^^ 

0.06(0.06) 


+0.43 


A 1 A 

-0.10 


A 1 A 
U. lU 






loo 
Is. 2 


07 


02 40 28.13(0.43) 


+61 15 50.8(2.6) 


0.07(0.04) 


0.12(0.09) 


+4.47 


+3.80 


1.62 


15.0 


14.5 


14.3 


08 


02 40 33.33(0.26) 


+61 18 05.8(1.9) 


0.12(0.04) 


0.24(0.10) 














09* 


02 40 33.74(0.20) 


+61 12 00.8(1.4) 


0.14(0.04) 


0.17(0.07) 


-0.36 


-0.70 


0.51 


19.2 


17.8 


16.7 


JQ*(a,b) 


02 40 44.94(0.05) 


+61 16 55.8(0.3) 


0.54(0.04) 


0.49(0.06) 


+0.07 


+0.20 


0.65 


6.4 


6.2 


6.2 


llM 


02 40 45.32(0.33) 


+61 12 42.3(2.0) 


0.10(0.04) 


0.17(0.09) 














12* 


02 40 52.15(0.06) 


+61 16 26.1(0.4) 


0.39(0.04) 


0.35(0.06) 


+0.43 


+0.10 


0.54 






18.1 


13*(a) 


02 41 07.76(0.06) 


+61 14 54.2(0.4) 


0.46(0.04) 


0.49(0.07) 


+0.79 


+0.30 


1.14 


20.5 


19.6 


18.2 


14* 


02 41 09.36(0.07) 


+61 14 34.5(0.5) 


0.34(0.04) 


0.34(0.07) 














15* 


02 41 09.94(0.03) 


+61 09 23.0(0.2) 


1.16(0.04) 


1.00(0.06) 














16* 


02 41 10.97(0.14) 


+61 12 58.1(1.2) 


0.18(0.04) 


0.23(0.08) 















(*) Sources already detected by M98. 

Candidate X-ray counterpart found for this source. 

Near infrared magnitudes taken from the 2MASS All-Sky Catalog of point sources. 



is detected in our Chandra observations. On the other hand and 
using the same Eq. # (1) criterion, VLA source # 11 is coincident 
with a Chandra X-ray source (labelled as # 27 in Tabled but no 
infrared counterpart is detected. If either the magnetar-like burst 
behaviour or the X-ray activity presenting QPOs are unrelated 
to LS I +61 303, then these radio sources appear to be potential 
alternative counterparts for these detections. 

We note that VLA source # 1 1 is located at the centre of the 
extended radio source D (a = 02''40"'45", 6 = +61°12:7) re- 
ported by M98 which was proposed to be a possible large-scale 
lobe powered by LS I +61 303. In our new map (Fig. [TJ, this 
object resembles more closely a background radio galaxy with 
bent bipolar jets. The X-ray detection might then originate in 
the central core of the radio galaxy. In this interpretation, the 
connection of this source to a magnetar-like flare/QPOs appears 
unlikely since this kind of behaviour is not typical of an extra- 
galactic object. 

The new VLA maps confirm all the extended radio emission 
features reported by M98. A possible supernova remnant origin 
was tentatively considered by these authors but the low surface 
brightness and absence of extended X-rays rendered it unlikely. 
Revisiting this scenario seems timely because the magnetar-like 
emitter, even if it is not LS I +61 303, is naturally expected to be 
a galactic compact object created during a past supernova event. 
Thus, we attempted to estimate the spectral index of the two 
brightest extended radio features located west {a - 02''40'"01', 
6 = +61° 13:9) and south (a = 02''40'"47'', 6 = +61°09:5) of 
LS I +61 303 in Fig. [T] (sources A and C in M98) by com- 
bining our 6 cm map with the 20 cm archive data also quoted 
in Table [1] In this process, the visibility data was appropriately 
constrained to match angular resolutions in both data sets. The 
resulting spectral index map is not conclusive enough due to the 
poor quality of the 20 cm data. Only the brightest parts of the 
western radio feature have a spectral index eiTor below 0. 1 sug- 
gesting a possible non-thermal origin (Sv v""^). However, 
Chandra observations reveal no X-ray counterpart, even for this 
western extended radio feature. Improved long wavelength ra- 
dio and X-ray observations are required to constrain more ac- 



curately the possibility of a faint supernova remnant associated 
with LS I +61 303 or any other source in the field. Concerning 
compact sources, only one source labelled as # 15 was detected 
at 20 cm. Its spectral index (-0.6 + 0.2) may suggest a non- 
thermal emission mechanism for this source. 

In conclusion, we have reported a handful of X-ray and radio 
sources within or close to the improved Swift-BAil location of 
the magnetar-like and QPO events towards LS I +61 303. No ob- 
ject displays any peculiar signature that could reveal a possible 
connection with these unusual phenomena. Although we cannot 
exclude them being alternative counterpart candidates, the most 
likely possibility is that LS I +61 303 is actually behind the ob- 
served flaring event. 
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Chandra 


Right Ascension 


Declination 


Rux 


Hardness 


Id. # 


(hms) 






(10-** phcm-2 s-') 


Ratio 


01 


02 39 54.899(0.024) 


+61 


12 39.80(0.24) 


0.66(0.27) 


0.82(0.08) 


02 


02 39 58.962(0.019) 


+61 


15 19.46(0.13) 


2.44(0.43) 


-0.93(0.35) 


03 


02 39 59.445(0.016) 


+61 


17 03.53(0.10) 


3.98(0.53) 


0.33(0.09) 


04 


02 40 01.019(0.024) 


+61 


16 45.73(0.14) 


2.43(0.43) 


-0.49(0.28) 


05 


02 40 02.592(0.023) 


+61 


17 28.46(0.16) 


2.69(0.48) 


0.21(0.15) 


06 


02 40 10.328(0.036) 


+61 


17 11.64(0.10) 


1.09(0.31) 


-0.48(0.46) 


07 


02 40 10.521(0.017) 


+61 


1144.04(0.13) 


0.77(0.27) 


-0.86(0.70) 


08 


02 40 10.804(0.014) 


+61 


12 16.98(0.13) 


1.01(0.29) 


-0.87(0.57) 


09 


02 40 16.604(0.018) 


+61 


17 37.79(0.18) 


1.22(0.34) 


-0.44(0.43) 


10 


02 40 22.166(0.030) 


+61 


17 57.30(0.16) 


1.09(0.33) 


-0.27(0.42) 


11 


02 40 22.802(0.018) 


+61 


08 47.61(0.26) 


0.83(0.29) 


-1.13(0.80) 


12 


02 40 22.821(0.007) 


+61 


15 30.59(0.05) 


2.42(0.40) 


0.02(0.17) 


13 


02 40 24.507(0.015) 


+61 


17 20.90(0.13) 


1.75(0.36) 


0.59(0.09) 


14 


02 40 26.729(0.011) 


+61 


16 20.20(0.13) 


1.37(0.33) 


-0.54(0.39) 


15 


02 40 26.859(0.016) 


+61 


16 29.57(0.15) 


1.36(0.33) 


0.28(0.19) 


16 


02 40 28.861(0.011) 


+61 


16 43.71(0.09) 


3.06(0.48) 


-0.44(0.24) 


17 


02 40 29.471(0.011) 
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